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ABSTRACT 



Aims. To model broad Ha wings observed in symbiotic binaries by an optically thin, bipolar stellar wind from their hot components as an 
alternative to that considering the Raman scattering of Ly/? photons on atomic hydrogen. 

Methods. Profile-fitting analysis. Comparison of the observed broad He wings and their luminosity with those predicted by the model. 
Results. Synthetic Ho- profiles fit excellently the observed wings for |Av| S 200km s -1 in our sample of 10 symbiotic stars during the quiescent 
as well as active phases. The wing profile formed in the stellar wind can be approximated by a function /(Av) oc Av~ 2 , which is of the same 
type as that arising from the Raman scattering. Therefore it is not possible to distinguish between these two processes only by modeling the 
line profile. Some observational characteristics of the Ha-emission, its relationship with the emission measure of the symbiotic nebula and a 
steep radio spectrum at 1.4 - 15 GHz suggest the ionized stellar wind from the hot component to be the dominant source contributing to the 
Ha wings during active phases. The model corresponding mass-loss rates from the hot components are of a few xl0~ 8 Mq yr~' and of a few 
x (10~ 7 - 10~ s ) Mq yr~' during quiescent and active phases, respectively. 
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• *"H 1. Introduction 



Ivan Winckel et al. I i 1993b and llvison et al~l (1 1994 presented 
a large survey of high- and low-resolution Ha line profiles 
of symbiotic stars. The profiles showed broad wings, in most 
cases extended well within the presented wavelength range of 
6540 - 6580 A. Their origin has been investigated by number 
of authors. A popular interpretation assumes Raman scattering 
of Ly/3 photons on atomic hydrogen to be responsible for filling 
in the broad Ha wings . This possibility was firstly pointed out 
bvlNussbaum er et ail i 19891) and the corr espo nding quantita- 
tive m odel was elaborated by iLeel (2000) and iLee & Hvung"! 
(2000). Other possibilities for the Ha wing formation mech- 
anism - rotating disks, electron scattering, fast stellar w ind 
and Ha damping wings - were also discussed. Robi nson et al. I 
(11994 modeled the Ha profiles on the assumption that they 
originate in an accretion disk. Acceptable fits were found only 
for CHCyg, AGDra and TCrB. Generally, the model wings 
were broader than the observed o nes. The possibility of the 
electr on scattering was analyzed bv lArrieta & Torres-P eimbert 
(2003) for a representative case of M2-9. They found unrealisti- 
cally high values of the electron temperature and concentration 
for the electron-scattering region. Concerning to the Ha damp- 
ing wi ngs there is no elaborated application for symbiotic bi- 
naries. Lee (2000) only discussed briefly this possibility for 



the case of SYMus (Sch mutz et al. lfl994 . He came to con- 
clusion that the wing emission arises in a much more ex tended 
region then that producing the line core. Skop al et al. I (|2002) 
modeled the extended Ha wings from active phases of CH Cyg 
by a spherically symmetric and optically thin stellar wind. A 
comparison between the modeled and observed profiles was 
satisfactory and also the derived mass-loss rate was in agree- 
ment with that suggested by the radio observations. Therefore 
we propose the fast stellar wind from the hot component in 
symbiotic binaries to be the most promising alternative to that 
considering the Raman scattering process. 

Accordingly, in Sect. 2 we introduce a model of a bipolar 
stellar wind at the optically thin limit to calculate the broad 
Ha wings. In Sect. 3 we compare our model profiles with those 
observed during quiescent and active phases of selected sym- 
biotic stars. In Sect. 4 we discuss observational characteristics 
of Ha profiles connected with the hot star wind. 

2. Model of the bipolar wind 

2.1. Signatures of the mass outflow 

Here we summarize main observational features of a mass 
outflow indicated for active symbiotic stars. They are: (i) 
Broadening of emission line profiles and/or the P-Cygni type 
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Fig. 1. Geometry of the wind model we used to calculate the line profile of Ha and its luminosity. The wind is produced by the 
central star (small circle at the mid of panels). A fraction of the wind is blocked by an optically thick disk/torus at the center that 
cut out a cone from the sphere with the opening angle n - 28q. By this way we simulate bipolar geometry of the stellar wind. The 
shadow belt on the side view (left) and/or the annulus from the frontal view (right) represents a part of the shell with the same 
radial velocity, — v(r) cos(8). Radiation from the annulus, cut out by 46 radians, is blocked by the central disk (the darker part). 
The angle 8 is counted in the plane containing the annulus. 



of profiles represent direct indications of a mass-flow from 
the star. Typical velocities are a few hundred of krns" 1 
(e.g. iFer nandez- Castro et al~lll995t Bjussbaum er et al.1 1 1 995t 
Skopal et al. 1997). (ii) A significant enhancement of the neb- 
ular emission in the continuum by a factor of ^10 relatively 
to quiescent phases ( Skopal 1 120051 Tables 3 and 4) can in 
part result from a supplement of new emitters into the particle 
bounded nebula, for example, due to an increase in the mass- 
loss rate from the active object, (iii) The radio light curves 
usually show a decline at beginnings of outbu rsts with an in- 
crease after the optical maxim um (e.g. Fernandez-Castro et al. I 
1 1 995t IBro cksopp et al. 2004). In some cases evolution from a 
point sourc e to a spatial structure was observed (e.g. C H Cyg 
and Z And: iKennv et al. Ill996l IBrocksopp etai~l 12004 . This 
can be a result of a high-velocity mass-outflow, which grad- 
ually becomes optically thin at radio wavelengths. Velocities 
from a few hundred of kms~' to 1000-^2000 km s" 1 can b e 
derived from images (e.g. CH Cyg, Crocker et al. l200lU2002l) . 
(iv) The X-ray emission can be also explained by interaction 
of the outflowing material with surrounding nebular gas. The 
extended X-ray emission in the CH C yg Chandra image was 
aligned with the optical and radio jets (Gal loway & Sokoloskil 
120041) . Thus the velocities connected with the X-ray emission 
can be similar to thos e derived from the radio. For ZAnd, 
ISokoloski et al. I (|2006) ascribed the X-ray emission from its 
major 2000-03 outburst to the shock-heated plasma as a conse- 
quence of the mass ejection from the white dwarf, (v) Emission 
line profiles of forbidden lines from highly ionized atoms can 



be produced by the wind-wind collisio n zone in a binary system 
(e.g. lWallerstein et al. Il984l Eriksson et al. 2004). Aspects of 
the wind-dynamics including collidin g winds in symbiotic bi- 
naries were reviewed bv lWalder & Foliriil lEoOO). 

Finally, we note that different features of the outflowing 
material in the spectrum can reflect different driving mecha- 
nisms. For example, ejected rings or shells produce broad pro- 
files with rather steep sides of all lines (FWHM ~ 0.7x FWZI). 
Classical novae 14 9 4 Aql and V475 Set d emonstrate this case 
dEvres et al.ll2005t IChochol et aTl 1200.4 . For a star with a 
spherically symmetric and optically thin wind, the line emis- 
sivity is proportional to the square of the particle concentration, 
which is diluted with the radial distance r as l/r 2 . In this case 
a strong line core, originating from the vicinity of the wind 
source, is accompanied with faint extended wings from large 
distances (FWHM «: FWZI). In the following model we con- 
sider this case. 

2.2. Geometry of the wind 

The geometry of our wind model is introduced in Fig. 1 . The 
model assumes an optically thin stellar wind with a spheri- 
cally symmetric str ucture a nd the origin at/around the central 
star. According to ISkopal I J2005I) we put an optically thick 
disk/torus, characterized with the height H and radius R D , to the 
center of the hot object. We assume the disk to be seen edge-on 
due to a high orbital inclination. The outer rim of the disk cut 
out the angle 26*o = 2 tan -1 (H/Ro) from a sphere with the cen- 
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ter at the accretor, and by this way simulates bipolar shape of 
the stellar wind with the opening angle n - 28q radians (Fig. 1 
left). The wind with this geometry produces line profiles that 
are symmetrical with respect to the reference wavelength (Ao). 
We assume that the particle concentration n(r) at any point in 
the wind is related to the mass loss rate M and the velocity v(r) 
via the mass continuity equation, i.e. 



n(r) = M/4nr /jmnv(r), 



(1) 



where fi is the mean molecular wei ght and mw is the m ass of the 
hydrogen atom. According to the Cast or et al. 1 ll 19751) model 
we approximate the velocity distribution in the hot star wind 
by 



v(r) = v m (\-RJr) 



(2) 



The velocity v(r) of the wind increases monotonically outward 
from its beginning at R w and asymptotically approaches the ter- 
minal speed Voo. The distance r is counted from the center of 
the star and the parameter ft characterizes an acceleration of 
the wind, i.e. the 'slope' o f v(r). A smaller /3 c orresponds to a 
faster transition to Voo (e.g. B ertout et al. Il985l) . 

2.3. Ha luminosity 

In our simplified approach we assume that the wind is fully ion- 
ized and completely optically thin in Ha with a steady mass- 
loss rate. The optically thin case is supported by the large ve- 
locity gradient in the stellar wind, because of its large termi- 
nal velocity (voo » v t h, where Vth is the thermal velocity). 
If a line photon, created by recombination in such a wind, 
has traveled a distance / > 2v t h/(dv/dP), it is Doppler shifted 
with respect to the surrounding gas by more than 2v t h and thus 
cannot be absorbed any mo re in the same line transition (e.g. 
lLamers. & Cassinellill999l) . Under such conditions the escape 
probability of the emitted photons will be close to 1 . However 
the optically thin condition can be attained only at large dis- 
tances from the source of the wind. A good agreement between 
the observed and modeled profiles for |Av] ^200kms _1 (see 
below, Fig. 3) suggests the validity of the optically thin regime 
from about 1.2-5-1.5 R w (Eq. (2) and parameters from Table 1). 

The total line luminosity, L(Ha), is related to the line emis- 
sivity of the wind, e a n e n + , by 



L(Ha) = 4ne, 



a J n e n 



(r)[l - w(r)] r 2 dr, 



(3) 



where r; is a certain distance from the source of the wind where 
the integration starts from (its quantity is given by the model; 
r\ > H, R w < H), s a = 3.56 x l(T 25 ergcm 3 s -1 is the volume 
emission coefficient in Ha for T e = 10 4 K. We assume it to be 
constant throughout the wind. n e and n + are concentrations of 
electrons and ions (protons). We assume a completely ionized 
medium (n e - « + ) and radial distribution of particles as given 
by Eq (1). The factor w(r) determines visibility of the wind for 
the outer observer. It can be expressed as 

w(r) = sinfl = H/r for r < ro, 
= sin6>o = H/ro for r > ro, 



(4) 



where r}. 



+ H 2 (Fig. 1). Therefore we have to integrate 
contributions from the shells between H and ro and those above 
the ro radius separately. Substitution of Eqs. (1), (2) and (4) 
into Eq. (3) and using dimensionless parameters x = R w /r, 
a = H/R w and / = R w /ro yields an expression for the Ha lu- 
minosity as 

M\ 2 



L(Ha) 



47r(//mH) : 



IMV- 1 

I Voo/ Rw 



X (/i + h), 



(5) 
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(7) 



x)V " (1-2/3) 
An alytical expression of the integral I2 was already introduced 
bv lSkopaletaLl J2002h . Thus comparing the Ha luminosity 
from observations to that predicted by Eq. (5) allow us to esti- 
mate the mass-loss rate, M. To determine theoretical values of 
L(Ha) requires a knowledge of parameters characterizing the 
wind: Voo, R v , j3 and 6q. To estimate their appropriate quantities 
we fit a synthetic-line profile to the observed one. This requires 
a different manner of integration of the wind's contributions. 
We introduce it in the following section. 

2.4. Ha profile 

To reconstruct the global line profile we redistribute the inde- 
pendent Doppler-shifted contributions from each volume el- 
ement of the wind expanding material in the radial velocity 
co-ordinates. The profile thus represents a 'broadening func- 
tion' resulting from a field of contributions which differ in 
emissivity and radial velocity. It can be compared only to 
high-velocity features in the profile produced by regions with 
r(Ha) < 1. A technique of integration can be understand with 
the help of Fig. 1 . Volume elements of the same radial velocity, 
RV - — v(r) cos 0, are represented by annuli around the line of 
sight, and can be expressed as 

A V = (2tt - 40)r 2 sin <9A6»Ar, (8) 

where the angle 48 corresponds to a fraction of the annuli, 
which radiation is blocked by the central disk. The 8 angle is 
counted in the plane containing the annulus of the radius r sin 8 
and is related to H and 8 as (Fig. 1 right) 

for H < r < ro, 



r sin6> 
sin 6>o 

sin# 



for r > ro. 



(9) 



A. Skopal: Broad Ha wings from the optically thin stellar wind of the hot components in symbiotic binaries 



6J0 
O 



9 H 
3 1 



o 
o 



RW Hya 
SY Mus 
CI Cyg 
YY Her 
Z And 
AX Per 
AG Peg 




-7.5 



-6.5 



log(dM/dt) [Ms^yr" 1 ] 



AR Pav 
AE Ara 
AX Per 
Z And 
PU Vul 



-5.5 



Fig. 2. Luminosity of Ha wings for |Av] > 200 km s , L a (2QQ), as a function of the mass-loss rate calculated according to 
Eq. (12) for parameters R w , f3 and Voo in Table 1. The observed quantities of Lq,(200) are denoted by filled diamonds (Table 1). 



In the sense of Eq. (3), radiative contributions of such cut-out 
annuli are AL(r, 0) = s a n(r) 2 AV. They correspond to a certain 
radial velocity Av (or AA) in the profile. With the aid of Eqs. (1) 
and (8) we can write them as 



(2tt -40) sin , 
AL(r, 0) = f- A0Ar [erg s" 1 ], 

where v(r) = v(r)/v„o and the factor 

f= g " (*)\ (11) 
(An iJ.mii) 2 V Voo ' 

Redistributing all the visible emissions according to their ra- 
dial velocities we obtain the resulting line profile. Integration 
of each shell begins from the direction at 6 = 6\ and ends at 
= n - 0i. Their contributions are summarized from H to a 
distance at which the wind's emission can be neglected. Thus 
the line luminosity can be expressed as 



H JBi 



' (27r-4f?)sin6l 
r 2 v(r) 2 



dftlr. 



(12) 



According to the relation (9), integration has to be divided into 
two parts: for H < r < ro, 0\ = sm~ l (H/r) and for r > ro, 
0\ = 8o. Contributions from the nearest regions to the wind 
origin, R w , are characterized with H < r <K ro and large val- 
ues of 0\. According to relations (1) and (2) they are strong 
due to a high density of the wind and contribute mainly to the 
line center, because of a small radial velocity. However, amount 
of their emission is reduced significantly by the large value of 
sin in our model (Eq. 9), i.e. a significant fraction of radia- 
tion is blocked by t he torus. I ntegration ends at a finite limit of 
200 R (Fig. CI in lSkopaletal. ll2002h . After summing of all 
the contributions we scaled the synthetic profile to the maxi- 
mum of the observed one. Examples are shown in Fig. 3. The 
input parameters, Vco and 0q can be inferred from observations, 
while the resulting parameters, R w , £ and f3, are given by an 
appropriate fit to the observed profile (Sect. 3.2). If we rewrite 
the relation (12) in the form 

L(Ua) = £ x h = And 2 F(Ha), (13) 



we can express the ratio 

F(Ha)\ 1 / 2 



— =8.5x 10 

Voo 



io ^(H«) j 



x d 



kms -1 



(14) 



(10) w here the distance d is in kpc and the observed bolometric flux 
F(Ha) in erg crrT 2 s~' . Following to Eq. (10) integral It, is given 
by the geometry of the wind and is proportional to the sum of 
all its visible volume elements. For wing fluxes, e.g. Fq.(200) 
(see Sect. 3.3), the 73 integral includes contributions only with 
the radial velocity |Av| > 200 km s~' . 



3. Comparison with observations 

3.1. Data sources 

In major part we used observations of Hg profiles from th e 
survey of Ivan Winckel etal.l dl993l) and llvison et al. I i 1994 . 
Other sources are referred in Table 2. If possible, we selected 
objects that were observed during both quiescence and activ- 
ity. We converted observed fluxes in relative units to fluxes in 
erg crrT 2 s^ 1 A -1 with the aid of simultaneous optical V and R 
photometry (Skopal et al. 1 120041 an d refe rences therein) and 
the model SED according tolskorjal (2005). Approximate cor- 
rections for emission lines dSkopal 1120031) were also included. 
T o deredden the fl uxes we used appropriate £b-v from Table 1 
ofSk^iD (2005). For the purpose of Sect. 4 we also estimated 
emission measure, EM^, of the symbiotic nebula at the dates 
of Ha observations (Table 2). For the sake of simplicity and 
availability we used dereddened fluxes, Fy, derived from U- 
magnitudes. Then according to Eq. (18) of Skopal (2005), 



EM obs ±And z — , 



(15) 



where sy is the volume emission coefficient per electron and 
per ion (erg cm 3 s~' A -1 ). We used the average value of s\j from 
both the sides of the Balmer jump corresponding to the electron 
temperature given by the SED ( Skop al 120051 . 
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Table 1. Parameters of Ha models (R w , ft, v x ), observed luminosities and corresponding mass-loss rates (M) and emission 
measures (EM W ) of the hot star wind. L „(200) denotes the luminosity of the Ha wings for |Av| > 200 kms -1 , while L a (0) is the 
total line luminosity. Distances are from SkopalN2005b . 



Object 


d 


Date 


R w 




Voo 


L a (0) 


L„(200) 


log(M) 






[kpc] 




[%] 




[kms" 1 ] 


[£©] 


[*©] 


[Moyr- 1 ] 


[cm" 3 ] 








Quiescent phases 










k w riya 


U.o 


1 C\l(Yl IQl 


0.040 


1.70 


1000 


i . i 


ft ^zL 


- / .Zj 


Z.U X 1U 


Z And 


1.5 


22/09/88 


0.040 


1.72 


1000 


5.1 


0.87 


-7.02 


5.4 xlO 58 


AX Per 


1.7 


01/11/93 


0.042 


1.70 


1000 


6.6 


1.2 


-6.96 


9.2 xlO 58 


CICyg 


2.0 


21/09/88 


0.042 


1.72 


1000 


8.0 


0.40 


-7.20 


3.2 xlO 58 






20/06/89 


0.040 


1.63 


1000 


3.3 


0.48 


-7.15 


2.6 xlO 58 


AG Peg 


0.8 


15/07/88 


0.042 


1.70 


1200 


4.6 


0.68 


-7.04 


4.4 xlO 58 


YYHer 


6.3 


20/06/89 


0.040 


1.75 


1000 


2.8 


0.40 


-7.20 


2.5 xlO 58 


SYMus 


1.0 


14/07/88 


0.042 


1.70 


1000 


3.0 


0.14 


-7.42 


1.1 xlO 58 










Active ph; 


ises 










ARPav 


4.9 


14/07/88 


0.90 


1.75 


1800 


39.5 


9.2 


-5.72 


5.2xl0 59 


AE Ara 


3.5 


14/07/88 


0.90 


1.80 


2000 


36.4 


6.3 


-5.78 


3.6xl0 59 


AX Per 


1.7 


22/09/88 


0.85 


1.80 


1600 


7.1 


2.0 


-6.07 


9.8xl0 58 


ZAnd 


1.5 


11/12/00 


1.70 


1.70 


2600 


20.8 


5.5 


-5.60 


2.9x1 59 


PUVul* 


3.2 


28/09/88 


1.67 


1.70 


2100 


144 


ll f 


-5.29 


1.5xl0 60 



d = 3.2 kpc, £ B _ V = 0.22lRud v"etal.ll999h . + = L„(500) 



3.2. Model parameters 

Geometrical parameters of our wind model are described in 
Sect. 2.2. Some limits for parameters Rd and H can be esti- 
mated from the effective radius, R 1 ^, of the hot star. This pa- 
rameter represents the radius of a sphere that produces the ob- 
served luminosity of the hot stellar source and can be derived 
from modeling the SED of the ultraviolet continuum dSkopal I 
120051) . In our model 4n(Rf) 2 = A nRnH. During quiescent 
phases we observe Rf * 0.15 Rq (Skopail Eooll Table 3). 
Assuming the ratio H/Ro = 0.1 then yields Rd = 0.5 Rq and 
H = 0.05 ^0. During active phases we assume flared disk with 
H/Ru = 0.3. Then parameters H and Ru are adjusted to the cor- 
responding R^ (a few of Rq) for objects we investigate here. 

The origin of the wind, R w , and ft in the wind law are 
model parameters, while v M is given by the extension of the 
wings. Values of /? w , and ft are critical for the synthetic pro- 
file. Generally, a larger value of ft corresponds to a slower and 
denser wind with a higher emissivity at a point r. Therefore the 
wind characterized with a larger value of ft produces a narrower 
profile. Regions close to the wind's origin (R w < H < r <K ro) 
have the largest emissivity, because of high densities for small 
values of both r and v(r); they contribute mainly to the line 
core. However, a fraction of their radiation is blocked by the 
outer rim of the disk in our model. In addition, optical prop- 
erties of these regions can deviate from the optically thin case. 
Therefore we do not aim to fit the core of the line by this proce- 
dure. At further distances (r 1.2 R w ) the wind is accelerated 
to v(r) ;> 200 - 250kms- 1 (Eq. (2), ft ~ 1.7). Contributions 
from these regions are more important to create the broad 
wings. A good agreement between the modeled and observed 
profiles for radial velocities |Av| S 200km s~' (Fig. 3) is con- 
sistent with our assumption on the optically thin regime from 
this distance above the source of the wind. 



3.3. Wing profiles, M and uncertainties 

We reconstructed synthetic profiles according to Eqs. (10) and 
(12). To obtain an appropriate solution we calculated a grid 
of models for reasonable ranges of R w and By this way 
we also estimated their uncertainties to 20% and 10%, respec- 
tively. A comparison of resulting models with observations is 
shown in Fig. 3 and the corresponding parameters are intro- 
duced in Table 1 . Models match perfectly the line profile for 
|Av| > 200kms~ 1 . Therefore to determine the mass-loss rate 
from the hot star, M, we compared the luminosity of wings for 
|Av| ^ 200km s~' (L„(200) in our notation) with that calcu- 
lated according to Eq. (14) for the same radial velocity inter- 
val. Dependences of L [V (200) on M for our models are plot- 
ted in Fig. 2. Uncertainties in the M values were determined 
from those of F (Ha), ft and R w . We estimated the uncertainty 
in the Ha fluxes to ~ 10%, which results mainly from deter- 
mination of the level of the local continuum (Sect. 3.1). As 
M oc ^LXHa), the uncertainty in fluxes increases that of M 
with only square root. We estimated them to 8% - 12%. 

4. Discussion 

We showed that the profiles from the ionized optically thin, 
bipolar stellar wind match well the observed Ha wings for 
I Av| ^ 200 km s _1 . The curve fitting the profile is proportional 
to Av~ 2 (Eqs. 10 and 12), which is of the same type as th at in- 
cludin g solely the Raman LyyS — > Ha scattering process jLeel 
2000). From this point of view it is not possible to distinguish 
contributions from the ionized wind and the Raman scattering 
in the wing profile directly. However, both the processes take 
place in very different regions of the binary. The ionized stel- 
lar wind in our model is located around the hot star, while the 
Raman-scattered photons originate in the neutral part of the 
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Fig. 3. Comparison of the modeled (M) and observed (O) Ha line profiles for selected symbiotic stars during their quiescent 
phases (top two raws of panels) and active phases (the bottom raw of panels). The systemic velocity was subtracted. The model is 
described in Sect. 2.4 and the corresponding parameters are in Table 1 . It fits the observed profiles for | Av| ^ 200 km s~' . Fluxes 
are in 10~ n ergcrrT 2 s~' A -1 . 



wind from the giant. Below we summarize observational char- 
acteristics of the Ha emission that could help to identify main 
sources of radiation contributing to its broad wings. 

4.1. Observed properties of Ha wings 

(i) During eclipses of the hot component by the giant the 
Ha emission in both the core and the wings decreased signifi- 
cantly. E xamples here are AX Per {skorjal etal.l200l[ Fig. 6) 
ARP av JOuiroga etal.ll2002l Fi g. 2), FN Sgr jBrandi etall 
l2005t Table 2. Fig. 4) and ZAnd Jskonal et al. 120061 Fig. 4). 



These observations suggest that a significant fraction of the 
broad Ha wings is formed nearby the hot star. 

(ii) Ilkeda&Tamural J2000h. lOuiroga et all i200& 
iMiko laievsk a et al. I d2003l) and Brandi et al. (2005} revealed 
that the radial velocities from the wings of Ha follow the 
orbital motion of the hot component in V1329Cyg, ARPav, 
AE Ara and FN Sgr, respectively. This implies that the region 
of the Ha wings formation is connected with the hot star. 

(iii) The wing profiles from our sample of objects are sym- 
metrically placed around Aq. In cases of CICyg (20/06/89, <p = 
0.85) and ARPav (17/07/88, <p = 0.71) we shifted the model by 
+ 12 and + 15 km s~', respectively, to match better observations. 
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Fig. 4. Total Ha line and wing luminosities, L a (Q) and ^(200), 
respectively, as a function of the observed emission measure 
EM t, s . The total line luminosities are divided by a factor of 5 
for a better visual display. Corresponding quantities are from 
Table 2. 



These shifts are consistent with the hot component orbital mo- 
tion. Only in the case of Z And (22/09/88, ip = 0.30) the model 
was shifted by + 10 km s , which is against the orbital motion. 
Nature of this difference is not clear. 

(iv) In our sample of Ha profiles we did not find systematic 
shifts, which could be associated with the Raman scattering. 
For example, the supposed blue-shifted component from the 
neutral wind at/aroun d the binary axis tha t moves against the 
incident Ly/3 photons llSchmid et al. Ill999|) and/or the redshift 
of the Ha wing center as suggested bv ljung & Lee 
were not indicated. From this point of view, a good agreement 
between systemic velocities determined independently from 
M-giant and Ha-wings radial velocities, respectively, for ob- 
jects referred in the point (ii), is consistent with our finding. 



4.2. Ha luminosity as a function of the nebular 
emission 

Here we investigate a relationship between the Ha luminosi- 
ties, L a (Q), Lq.(200), and the emission measure of the symbiotic 
nebula, Z?M bs> which is due to photoionization. Relevant data 
are described in Sect. 3.1, summarized in Table 2 and plotted 
in Fig. 4. 

During quiescent phases the nebular emission originates in 
the ionized part of the wind from t he giant as given by the STB 
llSeaquist. Tavlor & Button I Il984l) model. This w a s inde pen- 
dently supported by finding of Nussbaum er et al. I Jl988) that 
symbiotic objects fit well the CNO abundance ratios of normal 
red giants. Following analyses showed that th e model is ap- 
plicable for most of quiescent sy mbiotics Ce.g. lSeaquist et al.1 
1993; Mikolaiew2kaiIZ]E002). Within this model the ion- 
ized wind from the giant can contribute to only the Ha core 
emission, because of its small radial velocity dispersion. The 
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broad Ha wings thus have to be of a different nature. Figure 4 
plots the relevant quantities as a function of log(£'M b s ). The 
observed dependencies are in qualitative agreement with the 
model: The total Ha emission, L a (Q), is a strong function of 
EM b s , whereas the wing emission, L ff (200), shows only a faint 
dependence. Assuming that the wing emission originates in 
the hot stellar wind, the corresponding mass-loss rates are of 
a few xl0~ 8 Mq yr -1 (Table 1, Fig. 2). However, such the wind 
can produce only very small nebular emission. According to 
Eqs. (3) and (5) its emission measure, EM W = J v n e n + [1 - 
w(r)] dV, can be expressed as 



h+h 



h /My I 

h) 2 \ Voo / R w 



(16) 



Model parameters (Sect. 3.2 and Table 1) yield £M W 
10 58 cm 4 . Thus during quiescent phases 



£M W «k EM n] 



(17) 



because EM \, S as 10 59 cm 3 (Table 2). Qualitatively, l arger hot 
star luminosity gives rise to larger emission measure JSkopal I 
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120051 Table 3) and probably drives a stronger hot star wind. 
This would explain the L Q ,(200)/log(£'M o b s ) dependence. As 
EM W <k EM D ^, this dependence is faint and the wing emission 
cannot rival that from the ionized giant wind. Consequently, the 
radio emission satisfies that from the giant's wind (i.e. the STB 
model as mentioned above) without a detectable influence of 
the emission from the hot star wind. 

During active phases the wing emission increases with 
log(£'M bs) by a factor of ~ 4.4 faster than during quies- 
cence (Fig. 4). This supports the idea that the source of emis- 
sion producing the broad Ha wings is the ionized hydrogen. 
Our Ha-wing models then imply that its kinematics corre- 
sponds to a fast wind from the hot star at rates of a few 
x (1(T 7 - 10" 6 )M© yr -1 (Table 1). The wind produces a large 
emission measure, because it increases with M 2 (Eq. 16). Our 
models yield EM W ^ 10 59 cm~ 3 , which can rival the EM Q \, S 
(Tables 1 and 2). Thus during active phases 

£M W <, EM obs . (18) 

Such the large emission influences also the radio spectrum. 
For example, during the active phase of Z And the steepness 
of the radio continuum between 1.4 and 15 GHz was a factor 
of ~2 larger than during quiescence (Fig. 5), which implies 
optically thicker conditions. Therefore the cm radio emission 
from the activity could be attribute d to optically thick f - f 
emission from the hot star wind (cf. Seaau ist & Tavlorll9 92). 
However, we need more radio observations for active objects 
and to elaborate corresponding quantitative model taking into 
account their configuration, for example, that suggested by 
JSkonal 1200.4 Fig. 27). 

Finally, we note that the presence of a strong hot s tar wind 
during active phases is consistent with the finding of Skopal 
J2005h that the low-temperature nebula (LTN) in active sym- 
biotics has a high emissivity (N e ~ a few times 10 9 cirT 3 ) and 
is located around the hot star, because it is subject to eclipses. 
Therefore the LTN in active symbiotics may be attributed to the 
emission of the wind from the hot component. 

5. Concluding remarks 

In this paper we introduced a model of the optically thin, bipo- 
lar stellar wind from the hot components in symbiotic bina- 
ries. We derived an expression, which relates the Ha luminos- 
ity from the wind to the mass-loss rate (Eq. 5) and calculated 
the corresponding line profile with the aid of Eq. (10). We ap- 
plied the model to the observed Ha profiles for 10 symbiotic 
stars during their quiescent and active phases. Synthetic pro- 
files provide a good fit to wings for |Av| ^ 200km s _1 from the 
line center. According to Eq. (10) the wing profile can be ap- 
proximated by the curve /(Av) oc Av~ 2 , which is of the same 
type a s that resulting from the Raman scattering process jLeel 
2000). Therefore it is not possible to distinguish between con- 
tributions from the ionized hot stellar wind and that from the 
Raman scattering process by only modeling the line profile. To 
support the former possibility we investigated relationship be- 
tween the emission from the Ha wings and the emission mea- 
sure of the symbiotic nebula. 



Wavelength [cm] 
30 20 10 6 2 




9 9.5 10 

log(Frequency) [Hz] 



Fig. 5. Radio observations of Z And at cm wavelengths during 
quiescence and the recent activity. Steeper continuum spectrum 
during activity suggests optically thicker condition s than dur- 
ing quiescence. The data are from Brocksop p et al. 1 (12004). 



We found that during quiescent phases the dependence 
L ff (200)/log(£'M o bs) is faint (cf. Fig. 4). The wing emission is 
relatively very small (Eq. 17). If the wings originate from a 
fast hot star wind, the corresponding mass-loss rates are of a 
few xlO~ 8 M0yr~' (Table 1). However, it is difficult to indi- 
cate their emission independently by other observations. Even 
in the radio the wing emission has no detectable effect - the ra- 
dio spectra satisfy radiation of the ionized wind from the giant 
(the STB model). 

During active phases the L„(200)/log(£ , M o bs) relation is a 
factor of ~4.4 steeper than that from quiescence. In this case 
the wing emission can represent a significant fraction of the ob- 
served nebular emission (Eq. 18) and thus can affect the radio 
spectrum. At cm wavelengths the steep continuum spectrum 
during the recent activity of Z And (Fig. 5) is consistent with 
an optically thick / - / emission from the hot star wind. 

The L Q .(200)/log(£'M o b s ) relationship and other character- 
istics of the Ha profiles (Sect. 4.1) suggest that the ionized 
hydrogen located around the hot star in the form of a fast 
stellar wind is the dominant source of the emission in the 
±(200^-2000) km s broad Ha wings during active phases. 
The corresponding mass-loss rates are of a few x(10~ 7 - 
lO _6 )M0yr _1 (Table 1). This finding allow us to attribute the 
LTN emission in active symbiotics to that from the hot star 
wind. 

Acknowledgements. The author thanks the anonymous referee for in- 
spiring comments. This research was in part supported by a grant of 
the Slovak Academy of Sciences No. 2/4014/04. 

References 

Arrieta, A., & Torres-Peimbert, S. 2003, ApJS, 147, 97 
Bertout, C, Leitherer, C, Stahl, O., & Wolf, B. 1985, A&A, 
144, 87 



A. Skopal: Broad Ho- wings from the optically thin stellar wind of the hot components in symbiotic binaries 



9 



Brandi, E., Mikolajewska, J., Quiroga, C, et al. 2005, A&A, 
440, 239 

Brocksopp, C, Sokoloski, J. L., Kaiser, C, et al. 2004, 

MNRAS, 347, 430 
Castor, J. I., Abbott, D. C, & Klein R. I. 1975, ApJ, 195, 157 
Chochol, D., Katysheva, N. A., Pribulla, T., et al. 2005, Contrib. 

Astron. Obs. Skalnate Pleso, 35, 107 
Crocker, M. M, Davis R. J., Eyres S. P. S., et al. 2001, 

MNRAS, 326, 781 
Crocker, M. M., Davis R. J., Spencer, R. E., et al. 2002, 

MNRAS, 335, 1100 
Eriksson, M., Johansson, S., & Wahlgren, G. M. 2004, A&A, 

422, 987 

Eyres, S. P. S., Heywood, I., O'Brien, T. J., et al. 2005, 

MNRAS, 358, 1019 
Fernandez-Castro, T., Gonzalez-Riestra, R., Cassatella, A., 

Taylor, A. R., & Seaquist E. R. 1995, ApJ, 442, 366 
Galloway, D. K, & Sokoloski, J. L. 2004, ApJ, 613, L61 
Ikeda, Y., & Tamura, S. 2000, PASJ, 52, 589 
Ivison, R. J., Bode, M. E, & Meaburn, J. 1994, A&AS, 103, 

201 

Jung,Y-Ch., & Lee, H-W. 2004, MNRAS, 350, 580 

Kenny, H. T., Taylor, A. R., Frei, B. D., et al. 1996, in Radio 

Emission from the Stars and the Sun, ed. A. R. Taylor & J. 

M. Paredes, ASP Conf. Ser., 93, 197 
Lamers, H. J. G. L. M., & Cassinelli, L. P. 1999, Introduction 

to Stellar Winds, CUP, Cambridge 
Lee, H-W. 2000, ApJ, 541, L25 
Lee, H-W., & Hyung, S. 2000, ApJ, 530, L49 
Mikolajewska, J., Ivison, R. J., & Omont, A. 2002, Adv. Space 

Res., 30, 2045 

Mikolajewska, J., Quiroga, C, Brandi, E., et al. 2003, in 
Symbiotic Stars Probing Stellar Evolution, ed. R. L. M. 
Corradi, J. Mikolajewska & T. J. Mahoney, ASP Conf. Ser., 
303, 147 

Nussbaumer, H, Schild, H., Schmid, H. M., & Vogel, M. 1988, 
A&A, 198, 179 

Nussbaumer, H, Schmid, H. M., & Vogel, M. 1989, A&A, 211, 
L27 

Nussbaumer, H., Schmutz, W., & Vogel, M. 1995, A&A, 293, 
L13 

Quiroga, C, Mikolajewska, J., Brandi, E., et al. 2002, A&A, 
387, 139 

Robinson, K, Bode, M. E, Skopal, A., et al. 1994, MNRAS, 
269, 1 

Rudy, R. J., Meier, S. R., Rossano, G. S., et al. 1999, ApJS, 
121,533 

Schmid, H. M., Krautter, J., Appenzeller, I., et al. 1999, A&A, 
348, 950 

Schmutz, W., Schild, H, Miirset, U., & Schmid, H. M. 1994, 
A&A, 288,819 

Seaquist, E. R., Taylor, A. R., & Button, S. 1984, ApJ, 284, 202 
(STB) 

Seaquist, E. R., & Taylor, A. R. 1992, ApJ, 387, 624 
Seaquist, E. R., Krogulec, M., & Taylor, A. R. 1993, ApJ, 410, 
260 

Skopal, A. 2003, Baltic Astron. 12, 604 
Skopal, A. 2005, A&A, 440, 995 



Skopal, A., Vittone, A. A., Errico, L., et al. 1997, MNRAS, 
292, 703 

Skopal, A., Teodorani, M., Errico, L., et al. 2001, A&A, 367, 
199 

Skopal, A., Bode, M. E, Crocker, M. M., et al. 2002, MNRAS, 
335, 1109 

Skopal, A., Pribulla, T, Vanko, M., et. al. 2004, 
Contrib. Astron. Obs. Skalnate Pleso, 34, 45 [ArXiv: 
astro-ph/0402141| 

Skopal, A., Vittone, A. A., Errico, L., et al. 2006, A&A, 453, 
279 

Sokoloski, J. L., Kenyon, S. J., Espey, B. R., et al. 2006, ApJ, 
636, 1002 

Walder, R., & Folini, D. 2000, in Thermal and Ionization 
Aspects of Flows from Hot Stars: Observations and Theory, 
eds. H. J. G. L. M. Lamers, and A. Sapar, ASP Conf. Ser. 
204,331 

Wallerstein, G, Willson, L. A., Salzer, J., & Brugel, E. 1984, 
A&A, 133, 137 

van Winckel, H., Duerbeck, H. W., & Schwarz, H. E. 1993, 
A&AS, 102, 401 



